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Abstract 



The convergence of the perturbation expansion for the effective interaction 
to be used in shell-model calculations is investigated as function of the mass 
number A, from ^4 = 4 to ^ = 208. As the mass number increases, there 
are more intermediate states to sum over in each higher-order diagram which 
contributes to the effective interaction. Together with the fact that the en- 
ergy denominators in each diagram are smaller for larger mass numbers, these 
two effects could largely enhance higher-order contributions to the effective 
interaction, thereby deteriorating the order-by-order convergence of the ef- 
fective interaction. This effect is counterbalanced by the short range of the 
nucleon-nucleon interaction, which implies that its matrix elements are weaker 
for valence single-particle states in "large" nuclei with large mass number as 
compared to those in light nuclei. These effects are examined by comparing 
various mean values of the matrix elements. It turns out that the contri- 
butions from higher-order terms remain fairly stable as the mass number 
increases from A = 4 to A = 208. The implications for nuclear structure 

calculations are discussed. 

PACS number(s): 21.60.-n, 24.10.Cn 
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I. INTRODUCTION 



One of the long-standing problems in nuclear many-body theory has been the conver- 
gence of the perturbative expansion for the effective interaction H^q (or equally well that 
of effective operators) derived from realistic nuclear forces to be used in nuclear structure 
calculations. Conventionally, the various terms which appear in a perturbative expansion 
for i^gff are displayed by way of Feynman-Goldstone diagrams, examples of such diagrams 
are shown in Fig. |T[ It is well known that a realistic nucleon-nucleon interaction V^-^ con- 
tains strong components of short range, which renders a perturbative expansion in terms 
of Vj^-^ meaningless. To overcome this problem, one takes into account the short-range 
correlations through the solution of the Bethe-Brueckner-Goldstone equation and considers 
a perturbation in terms of the nuclear reaction matrix G. The wavy lines in Fig. |l| rep- 
resent such G interactions. However, higher-order perturbative contributions in terms of 
the G-matrix, may be large, and the convergence of the perturbative expansion slow or not 
convergent at all. Actually, Barrett and Kirson showed that third-order contributions to 
the effective interaction were substantial, and raised the question whether the perturbative 
Rayleigh-Schrodinger (RS) expansion in terms of the nuclear G-matrix did converge at all. 
Schucan and Weidenmiiller even pointed out that the order-by-order expansion of the ef- 
fective interaction will ultimately diverge when so-called intruder states are present. Typical 
intruder states for nuclei like ^^O and ^^Ca are four-particle- two- hole core-deformed states. 
It ought however to be mentioned that for nuclei with more valence nucleons in e.g. the 
oxygen mass area, such intruder state configurations may not be important, and a two-body 
(or many-body) effective interaction defined within the lsO(i-shell only, may represent the 
relevant degrees of freedom. 

Most microscopic investigations of H^q have been performed for nuclei in the IsOd- 
shell, with few valence nucleons outside a ^^O core. However, when one extends the area of 
investigation to nuclei in the mass regions of calcium, tin or lead, one has to face the problem 
that for diagrams like those displayed in Fig. ^, there are more intermediate states which 
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contribute to each diagram of the effective interaction in e.g. the IpO /-shell than in the 
IsOd-shell. Moreover, the energy spacing between the various major shells is also smaller for 
nuclei in the IpO /-shell than for those around ^^O. This leads to smaller energy denominators 
which should enhance third-order or higher-order contributions. Thus, the combined action 
of the above effects could seriously deteriorate the order-by-order convergence (if it does 
converge) of the effective interaction. The only mechanism which could quench these effects, 
is the fact that the matrix elements of G calculated in the lpO/-shell should in general 
be weaker than those in the IsOrf-shell. The single-particle wave functions for the states 
around the Fermi energy exhibit larger radii and, as the nucleon-nucleon interaction is of 
short range, the matrix elements of G should be weaker for the heavier nuclei. The same 
arguments apply of course as well for the tin and lead regions. 

It is then the scope of this work to study the convergence of the effective interaction 
in terms of the mass number A, in order to assess whether higher-order contributions to 
the two-body effective interaction decrease or increase as A increases. To achieve this, we 
calculate all non-folded valence linked diagrams through third-order in the interaction G, 
and sum higher-order folded diagrams to infinite order for the mass regions beyond closed- 
shell cores with A = 4, A = 16, A = 40, A = 100 and A = 208. The details on how to 
obtain these effective interactions are briefly sketched in the next section, together with our 
results and discussions. Some concluding remarks are given in section three. 

II. COMPUTATIONAL DETAILS AND RESULTS 

There are basically two main approaches in perturbation theory used to define an effec- 
tive operator and effective interaction, each with its hierarchy of sub-approaches. One of 
these main approaches is an energy-dependent approach, known as Brillouin-Wigner per- 
turbation theory, while the Rayleigh-Schrodinger (RS) perturbation expansion stands for 
the energy independent approach. The latter is the most commonly used approach in the 
literature an approach which we will also employ here. It is then common practice in 
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perturbation theory to reduce the infinitely many degrees of freedom of the Hilbert space 
to those represented by a physically motivated subspace, the shell-model valence space. In 
such truncations of the Hilbert space, the notions of a projection operator P on the model 
space and its complement Q are introduced. The projection operators defining the model 
and excluded spaces are defined by 

i=l 

and 

oo 

Q= E mi^il, (2) 

i=D+l 

with D being the dimension of the model space, and PQ — 0, P'^ — P, Q'^ — Q and 
P + Q = I. The wave functions are eigenfunctions of the unperturbed hamiltonian 
Hq — T + U (with eigenvalues £i), where T is the kinetic energy and U an appropriately 
chosen one-body potential, in this work that of the harmonic oscillator (h.o.). The oscillator 
energies hil will be derived from hil = 45A~^/^ — 25A~^/^, A being the mass number. This 
yields HQ = 18.4, hfl = 13.9, HQ = 11.0, HQ = 8.5 and HQ = 6.9 MeV for A = 4, A = 16, 
A — 40, A — 100 and A — 208, respectively. The full hamiltonian is then rewritten as 
H = Hq + Hi with Hi = V^-^ — U, Vj^-^ being the nucleon-nucleon (NN) interaction. 
Below we will replace V^-^ by the G-matrix, which will be used as the starting point for 
our perturbative treatment. 

Following the above philosophy, we choose the model spaces which are believed, from 
both experiment and theoretical calculations, to be relevant for calculations of particle- 
particle effective interactions in the mass areas from A — A to A — 208. These are the 
0^3/2 and Opi/2 orbits for -f^gg in the mass area of A = A, the 0^5/2, 0^3/2 and lsi/2 orbits 
for A — 16, the 1^3/2, lpi/2, O/7/2 and O/5/2 orbits for nuclei in the mass region of A = 40 
and the O/111/2, 0^7/2, lc?5/2, lds/2 and 2si/2 orbits for A — 100. For these systems, the 
closed-shell cores ("^He, ^^O, "^"Ca and ^''°Sn) have equal numbers of protons and neutrons, 
and the model spaces are the same for both protons and neutrons. For lead however, with 
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Z = 82 and = 126, the proton and neutron model spaces are different, i.e. the orbits 
0«i3/2, O/19/2, 1/7/2, 1/5/2, 2p3/2 and 2pi/2 for the proton model space and Ozn/a, 0ji5/2, lfi^9/2, 
1(^7/2, 2(^5/2, 2(i3/2 and 3si/2 for the neutron model space. Since the effective interaction 
theory we will employ is tailored to degenerate model spaces, we will make no attempt to 
derive for lead an effective proton-neutron interaction for these two model spaces. Moreover, 
as discussed in Ref . , a multishell effective interaction may show strong non-hermiticities, 
or even divergencies if a h.o. basis is used. Thus, for A = 4 to A = 100 we will discuss both 
isospin r = and T = 1 effective interactions, whereas for lead we restrict the attention to 
Tz = —1 and Tz = 1, where Tz is the projection of the total isospin. 

For the above model spaces, there are in total 15 matrix elements for the effective in- 
teraction of A = 4, 63 for A = 16, 195 for A = 40, 353 for A = 100, 711 for the neutron 
model space of A = 208 and 353 for the proton model space of A = 208. The effective 
interactions for A = 16, ^4 = 40 and A = 100 are listed in Ref. 0, and have been tested in 
nuclear structure calculations and a good agreement with the experimental data obtained 
for several isotopes in these mass areas. The spectra for isotopes in the lead region will be 
published elsewhere 

Having defined the various model spaces, the next step in our calculation is to obtain 
the nuclear reaction matrix G, given by 

UJ — Uq 

where uj is the unperturbed energy of the interacting nucleons, and Hq is the unperturbed 
hamiltonian. For the bare NN interaction we use the One-Boson-Exchange potential Bonn 
A defined in Table A.l of Ref. The operator Q is a projection operator which prevents 
the interacting nucleons from scattering into states occupied by other nucleons. Note that 
the exclusion operator used in the calculation of the G-matrix in this work is different from 
the Q operator used in the evaluation of the effective interaction. The definition of the Pauli 
operator for the G-matrix can be found in Refs. where the so-called double-partitioned 
scheme has been used. This means that low-lying two-particle states are excluded by Q from 

6 



the intermediate states in the Bethe-Goldstone Eq. (|^). For the example of the IsOcZ-shell 
this exclusion refers to states with two nucleons in the IpO /-shell. As a consequence, we 
have to include in our perturbation expansion ladder type diagrams, such as (2-3) in Fig. |1|, 
where the allowed intermediate states are those of the IpO /-shell or corresponding ones for 
the other model-spaces. 

The next step is to define the so-called Q-boxS given by 

PQP = PHiP + P (hi—^Hi + Hi—^Hi—^Hi + ..)P, (4) 

where we will replace H\ with G [G replaces the free NN interaction Vj>^]v^). The Q-box 
is made up of non-folded diagrams which are irreducible and valence linked. A diagram 
is said to be irreducible if between each pair of vertices there is at least one hole state or 
a particle state outside the model space. In a valence-linked diagram the interactions are 
linked (via fermion lines) to at least one valence line. Note that a valence- linked diagram 
can be either connected (consisting of a single piece) or disconnected. In the final expansion 
including folded diagrams as well, the disconnected diagrams are found to cancel out 0. 
This corresponds to the cancellation of unlinked diagrams in the Goldstone expansion 
We can then obtain an effective interaction H^i^ = Hq + V^q in terms of the Q-box, with 
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Observe also that the effective interaction V^^ is evaluated at a given model space energy 
uj, as is the case for the G-matrix as well. For all mass areas, we fix a; = —20 MeV. The 
first iteration is then given by 

^eff = <3- (6) 

In this work we define the Q-hox to consist of all diagrams through third order in the G- 
matrix, as discussed in Ref. |^. Less than ten iterations were needed in order to obtain a 



^Not to be confused with the Pauli operators Q of the effective interaction and the G-matrix Q. 
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converged effective interaction for the various values of A. For further details, see Ref. 0. 
In the calculation of the various diagrams, we limit the intermediate state excitations to 2hu: 
in oscillator energy, an approximation which is viable if one employs an NN potential with 
a weak tensor force (such as the Bonn potential used here), as discussed by Sommermann et 
al. 0. It is the aim of this study to explore the effects of the various contributions to V^g. 
As it will be rather confusing to discuss the effects for individual matrix elements (recall 
that depending on the model-space there are up to few hundred matrix elements), we define 
averages of matrix elements by 

5:5:(2j+i)(2r+i) {j\o\j)j^ 

^^^'"^ ^ ~ ^^(2J+l)(2T+l) ' 

j JT 

where the summation index j refers to all two-particle states of the model-space under 
consideration, coupled to angular momentum J and isospin In the averaging procedure 
defined in this equation we have weighted the matrix elements by the factor (2J -|- 1)(2T -|- 
1) since this factor accounts for the degeneracy of two-particle states with respect to the 
projection quantum numbers and occurs e.g. in the calculation of the energy if all valence 
states are occupied. It turned out, however, that the main features of the results discussed 
below are obtained as well, if this weighting factor is dropped. For the operator O we will 
consider V^l), which corresponds to the bare G matrix, Q^, the Q-hax including terms 
up to second order in G without folded diagrams, and V^'^'^ (V^*-^-*) the effective interaction 
including all Q-box diagrams up to second (third) order plus all folded diagrams derived 
from these Q-boxes. Note, that the average defined in Eq. (|^) includes only diagonal matrix 
elements. In order to study if the conclusions remain valid for all matrix elements we also 
define a mean value including all matrix elements by 



^In tables | and || we omit to divide with the number of configurations, as this gives rather small 
numbers for the heavier nuclei. 
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EE(2^ + l)(2T + l)(^|O|0jT 

" E5](2j + i)(2r+i) ' 

fci JT 

where the summation indices k and / include again all two-particle states of the model- 
space considered. Beside these averages, which include matrix elements of isospin T = 
and T = 1, we will also report on results where the averaging is restricted to one of these 
isospins only. 

Results for the mean values of diagonal matrix elements (see Eq. (|^)) are listed in table |, 
while averages including the non-diagonal matrix elements as well (see Eq. (^) are presented 
in table Jl| for the various model-spaces considered. 

Inspecting these tables one observes very clearly that the mean values for the matrix 
elements are getting less attractive for the model spaces referring to heavy nuclei. This trend 
can be observed independent on the approximation used to calculate l^g. This behavior 
reflects the fact that also the effective interaction, calculated with inclusion of higher order 
terms, is of short range and therefore, as we discussed already above, yield weaker matrix 
elements for the valence nucleons in heavy nuclei as compared to the light systems. 

Furthermore, we observe some features which are valid independent on the mass number 
and model space considered: 

• The inclusion of second-order Q-hox diagrams in Q*^^) yields a substantial attraction 
for the T = matrix elements and a repulsion for T = 1. This difference may be 
understood by the following argument: For the T = 1 channel, the major mechanism 
which accounts for the difference between first and second order, is provided by the 
core-polarization diagram in (2-2) of Fig. |l]. Moreover, in the T = 1 channel, the tensor 
force component of the nucleon-nucleon interaction is not so important, whereas in the 
T = channel the contribution from the '^Si-'^Di partial wave plays an important role 
in ladder-type diagrams, such as several of the folded diagrams, or the particle-particle 
ladder diagram in (2-3) of Fig. |I]. Typically, for many J = 1 and T = particle-particle 
effective interactions, the particle-particle ladder is of the size of or larger than the 
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core-polarization diagram, while for J = and T = 1, the core-polarization diagram 
and the G-matrix yield the largest contribution to the effective interaction. 

• The inclusion of folded diagrams yields a repulsive trend going 

formQ(2) to 1/(2) _ rpj^g 

effect is again much larger in the T = than in the T = 1 matrix elements, which can 
as well be understood from the importance of the particle-particle ladder diagrams in 
the T = states. Comparing the results of V^^^ and V^'^^ one observes a repulsion for 
both isospins. 

• Contrary to this repulsion due to the second-order terms in the folded-diagram expan- 
sion, the additional inclusion of terms of third order in G yields some attraction in 
V^^^ as compared to V^'^^ . Except for the case of ^He, the effect of third-order terms 
is very weak for the T = 1 states. This was also observed in Ref. in the study of 
the spectra of nuclei with valence particles being only neutrons or protons. There the 
authors noted that the spectra of e.g. ^'^O or ^^Ca obtained with either a second-order 
or third-order effective interaction were quite similar. For calculations of the effective 
interaction for lead or tin, this is a gratifying property since it means that one needs 
only to evaluate the Q-box to second order and sum all folded diagrams. 

Finally, in order to discuss the convergence of the perturbation expansion, we compare 
in table |T| the ratios evaluated from the mean values defined in Eq. These ratios reflect 
of course the same features which we already discussed above. They emphasize, however, 
in a much better way that the different ratios are rather insensitive on the mass number 
which is considered. This means that one can expect the convergence of the perturbation 
expansion for the residual interaction to be as good (or bad) for heavy nuclei as for the light 
nuclei around ^^O, which are usually studied. 

For nuclear structure studies of heavy nuclei with neutron numbers quite different from 
the proton number one typically considers model-spaces, which are separate for protons and 
neutrons, ignoring the residual interaction beyond the mean-field approximation. For these 
cases (isospin T = 1), the effects of terms of second order in G seem to be rather important 
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with y*^^) containing a correction of around 50 percent of the average of V^^\ However, it 
is encouraging to note that the inclusion of third order terms yields a correction of only 5 
percent or even below. 



III. CONCLUSIONS 

We have studied the behavior of the perturbation expansion for the effective interaction 
to be used in shell-model studies of nuclei with various mass numbers. Inspecting appropriate 
mean values of matrix elements, we have found that the fact that the G-matrix becomes 
smaller in absolute value with increasing mass numbers, counterbalances the effects that 
there are more intermediate states to sum over and that the energy denominators become 
smaller in each individual diagram of the effective interaction. Therefore, the convergence 
of the perturbation expansion seems to be rather insensitive to the nuclear mass number. 
We observe that various features of the folded-diagram expansion, which had been discussed 
for the mass region A ^ 16, can also be found in heavy nuclei. The nuclear structure 
calculations for heavy nuclei are mainly sensitive to the proton-proton and neutron-neutron 
residual interactions. For these T = 1 matrix elements the third-order and second-order 
averages are very close, indicating that for this isospin channel one can approximate the 
effective interaction by including all diagrams to second order plus folded diagrams to all 
orders. For T = 0, one still needs to account for third-order contributions. 

The fact that third-order contributions seem to stabilize for heavier nuclei, has also 
important consequences for nuclear structure calculations in nuclei in the mass regions of 
e.g. ^'^^Sn and ^'^'^Pb. This means that the methods used to calculate the effective interaction 
for valence nucleons, applied mainly in the mass regions of ^^O and ^°Ca, can be applied to 
the mass regions of ^^^Sn and ^°^Pb, as done recently in Refs. ||6|JT0[ . 
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TABLES 

TABLE I. The mean values for diagonal matrix elements calculated according to Eq. (7) in 
model-spaces with cores as indicated in the first row assuming various approximations for the 
effective interaction. Averages are listed for all isospins {J2t) ^ '^^11 as for T = {J2t=o) 
T = 1 {J2t=i)- I^or lead, results for averages in the proton-proton model-space {J2pp) ^^id the 
neutron-neutron model-space {J^nn) listed. All entries in MeV. 



J2t=o 



^He 

-3.42 -5.98 -2.23 

(QP)) .3.38 -7.04 -1.68 

.2.87 -5.85 -1.48 

.3.42 .7.21 -1.65 

160 

-1.39 -2.67 -0.88 

(Q(2)) .1.24 -3.10 -0.50 

.0.99 .2.39 -0.43 

(1^(3)^ _2.80 -0.44 

40Ca 

-0.68 -1.35 -0.43 

.0.58 -1.57 -0.21 

.0.45 .1.20 -0.18 

.0.53 -1.51 -0.17 



100 



Sn 
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(F(3)) 

208 p]^ 

(yd)) 

(F(2)) 
(F(3)) 



-0.28 
-0.23 
-0.19 
-0.21 



-0.57 
-0.61 
-0.48 
-0.57 

-0.09 
-0.04 
-0.04 
-0.04 



-0.17 
-0.10 
-0.08 
-0.07 

-0.06 
-0.04 
-0.03 
-0.03 
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TABLE II. Mean values for diagonal and non-diagonal matrix elements calculated according 
to Eq. (8). Further details see caption of table I 



El 



Y,T=0 



^He 



-2.56 
-2.59 
-2.21 
-2.61 



-3.06 
-3.56 
-2.95 
-3.69 



-2.32 
-2.13 
-1.84 
-2.09 



160 



-0.72 
-0.68 
-0.54 
-0.59 



-1.13 
-1.29 
-0.99 
-1.16 



-0.56 
-0.44 
-0.36 
-0.37 



40Ca 

(F(3)) 



-0.28 
-0.27 
-0.21 
-0.23 



-0.41 
-0.48 
-0.36 
-0.46 



-0.22 
-0.19 
-0.15 
-0.14 



-0.13 
-0.12 
-0.09 
-0.11 



-0.30 
-0.33 
-0.26 
-0.30 



-0.07 
-0.04 
-0.04 
-0.04 
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208pb 

(yd)) 



-0.03 
-0.01 
-0.01 
-0.01 



-0.02 
-0.01 
-0.01 
-0.01 



17 



TABLE III. Ratios of mean values for diagonal matrix elements calculated according to Eq. 
(7). These ratios have been evaluated from mean values with better precision than those listed in 
table I. 



^He 

Er 0.84 1.19 

Er=o 0-98 1.23 

Er=i 0-66 1.11 

160 

Et 0-71 1.13 

Et=o 0-90 1-17 

Er=i 0-49 1.02 

40Ca 

Er 0.66 1.18 

Et=o 0.89 1.26 

Er=i 0.42 0.94 

looSn 

Et 0.68 1.11 

Er=o 0.85 1.19 

Er=i 0.49 0.94 

208 p]^ 

Epp 0.44 1.00 
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FIGURES 



FIG. 1. Different types of valence-linked diagrams. Diagram (2-1) is the G-matrix, diagrams 
(2-2)-(2-4) are second-order terms in G, while diagrams (2-5)-(2-8) are examples of third-order 
diagrams. 



20 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9509006vl 



